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ABSTRACT 

We have done a statistical analysis of Very Long Baseline Array (VLBA) data of water masers 
in the star-forming regions (SFRs) Cepheus A and W75 N, using correlation functions to study the 
spatial clustering and Doppler-velocity distribution of these masers. Two-point spatial correlation 
functions show a characteristic scale size for clusters of water maser spots <1 AU, similar to the 
values found in other SFRs. This suggests that the scale for water maser excitation tends to be <1 
AU. Velocity correlation functions show power-law dependences with indices that can be explained 
by regular velocity fields, such as expansion and/or rotation. These velocity fields are similar to 
those indicated by the water maser proper-motion measurements; therefore, the velocity correlation 
functions appear to reveal the organized motion of water maser spots on scales larger than 1 AU. 
Subject headings: ISM : individual (Cepheus A) — ISM : individual (W75 N) — ISM : jets and outflows 
— masers — methods : statistical 
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1. INTRODUCTION 

Water maser emission is commonly found in star- 
forming regions (SFRs) within clumps of molecular gas 
with high densities 10^ — 10^ c m~'^) and warm tem- 
peratures (~400 K) (jElitzur et al.ill989D . These physical 
conditions can be reached in both the inner parts of ac- 
cretion disks around young stellar objects (YSOs) and 
the gas compressed by shock waves generated by out- 
flows associated with these objects. 

Water masers are very compact and intense (~ 
10^'^ cm, Tb ~ 10"'^^ K), allowing studies through Very 
Long Baseline Interferometry (VLBI) observations to an- 
alyze in detail the spatio-kinematical distribution of the 
gas around YSOs with very high angular and spectral 
resolutions. In particular, sources with a large number 
of masers allow us to obtain through a statistical analy- 
sis the clustering and velocity distribution of masers from 
scales of a few hundred of AUs down to less than 1 AU. 
Statistical studies give important information about the 
typical size for clusters of masers that may be related to 
a fundamental scale of its excitation, as well as the sta- 
tistical properties of both the clustering and the velocity 
field. 

In this context. iGwinn et al. I (1992) carried out VLBI 
water maser observations of the SFR W49. They identi- 
fied 271 maser features with 625 maser spots, which were 
used for a statistical study. A maser spot is an individ- 
ual component of the maser emission occurring at a given 
velocity channel and position. Maser spots tend to clus- 
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ter in positi on and Doppler velocity, typically 1 AU and 
0.5 km s-i (G winn "li^Ma"). Such clusters or groups of 
maser spots are called maser features. Physically, maser 
features should be small clouds sup porting populat ion in- 
version by a pumping mechanism (|Gwinn 1 ll994bl) . The 
statistical analysis of spatial and velocity distributions 
for both spots and features reveals that turbul ent mo- 
tions dominate on a spatial scale of ~ 1-300 AU (|G wmn I 
19944 

ilmai et al. I (|2002l) carried out a similar study using 
Very Long Baseline Array (VLBA) maser data of the 
SFR W3 IRS 5. For the statistics, Ilmai et al. I pOOl 
used 905 maser spots grouped in 152 maser features. 
They found that spots form features with a typical size 
of ~0.5 AU. The statistical analysis of the Doppler ve- 
locities, specifically, the velocity correlation functions for 
maser spots follow a power-law dependence in the range 
of 0.04-300 AU with an index of ~0.29, that is con- 
sistent with the Kolmogorov value of 1/3, the value ex- 
pect edfor^n^om^ressiMeflA^ a turbulent velocity 
field (|Kolmogorovl 119411 : iStrelnit ski 2007 ). Simi l ar re - 
sults have been also obtained by IStrelnitski et al.l (|2002[ ) 
in other five SFRs. 

VLBA water maser observations toward the SFRs 
Cepheus A and W75 N revealed rem arkable mi- 
crostructures (jTorrelles et al. |[200TH I2003D . These mi- 
crostructures exhibit a coherent an d well-ordered s patio- 
kinematical behavior at AU scales (jUscanga et al.i r2005') . 
Proper-motion measurements of water masers suggest 
the presence of organized motions of stru ctures with sizes 
from tens to a few hundreds of AUs (jTorrelles et al. I 
l2001bl 120031 ). Here we study the spatial and velocity 
distribution of the water masers in these two SFRs using 
a statistical analysis to investigate whether organized or 
turbulent motions dominate over spatial scales from a 
few hundred of AUs down to less than 1 AU. 

This paper is organized as follows. In Section 2, we 
define the correlation functions used in this statistical 
study. In Sections 3 and 4, we give a brief introduction to 
the SFRs, Cepheus A and W75 N, and present the results 
of the correlation functions, respectively. In Section 5, 
we discuss the clustering of water maser spots and give a 
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possible interpretation for the velocity field traced by the 
water masers, based on the statistical analysis. Finally, 
we summarize our conclusions in Section 6. 

2. STATISTICAL ANALYSIS OF MASER SPOTS 
2.1. Two-point spatial correlation function 
We ado pt the defini t ion p reviously used by iWalker I 



([T98¥) and llmai et al. I ()2002[ ). considering the two-point 
spatial correlation function of maser spots as the number 
of spots per unit angular area with a given separation 
on the plane of sky, Ar, from an arbitrary spot. This 
function can be expressed as 



ns{Ar)dn 



^spot 



(1) 



where rigpot is the total number of maser spots, the vec- 
tors Ti and r j determine the positions of the spots, and 
the indices i, j run over all spots. Here 



ns{r) 



1 when Ar < r < Ar 
otherwise , 



dr . 



(2) 



where Ar = {Ax^ + Ay'^y/^, Ax and Ay are the off- 
sets in the right ascension and declination axis direc- 
tions, respectively; dr is a value for the separation of 
the successive bins in the rig vs. Ar plot, and (if2 = 
2TrAr dr + 7r((ir)^ is the total area of the annulus where 
the counting is done. 

Frequently, the two-point correlation functions can be 
approximated by a power law, ns{Ar) — noAr", on cer- 
tain ranges. A uniform o r random distribution would 
produce an index a = 0.0 () Walker 1 119841 : iGwinn 1 119945 
llmai et al. l[200l . 

The two-point spatial correlation function allows us to 
find a characteristic scale size for clusters of maser spots 
as well as the statistical properties of the clustering of 
spots. 

2.2. Velocity correlation functions 

Velocity correlation functions measure the variation in 
Doppler velocity with spot separation. We have used th e 
velocity correlation functions defi ned bvlWalkerl ()1984D . 
iGwinn I Hl994aD . and llmai et al. I ()2002t) in their statisti- 
cal studies. 

The first one is the rms Doppler-velocity difference as 
a function of spot separation (Ar) on the plane of sky. 
It is expressed by 



^.(Ar) 



E.A^^-VJ)'M\r^-'r,\) 
J2ijns{\ri-rj\) 



1/2 



(3) 



where ns{r) is given by equation (2). The second one is 
the median Doppler-velocity difference of the spots as a 
function of Ar. It can be written as 



M^(Ar) = Med(|Fj - V, |) for i, j 

so that Ar < |r,; — Tj | < Ar 



dr . 



(4) 



Again, the indices i, j run over all spots, Ar = {Ax^ + 
Ay'^y/^, and dr is the separation of the adjacent bins in 
the Vs and Mg vs. Ar plots. Generally, these velocity 
correlation functions can be fitted by a power law (Ar") 
in the whole range. The Ms function is less sensitive to 



outlying points with a large velocity difference ()Gwinn I 
Il994a). 

Velocity correlation functions give information about 
the gas kinematics traced by maser spots. 

3. STAR-FORMING REGIONS: CEPHEUS A AND W75 N 

Cepheus A is a high mass SFR located at a distance 
of ~700 pc (|Moscadelh et al.ll2009[ ). The brightest radio 
cont inuum source of the region, HW2, is a thermal radio 
let fHughes & Wouterloofl [l9M; iRodriguez e t al.l 119941 : 
[Gomoz et al. 1999). This source is associated with aclus- 
ter of water masers detected by iTorrelles etaP (|1996l ) 
using the Very Large Array (VLA). Through VLB A 
multi-epo ch water maser observa tions toward the radio 
jet HW2. ITorrelles et al. I (|2001bD detected ~1000 water 
maser spots in each of the three observed epochs, form- 
ing remarkable linear/arcuate structures with sizes of the 
order of tens of AU. These structures are grouped in five 
subregions (named as Rl to R5) within --0'.'6 (~400 AU) 
from HW2. The most interesting structure is an arc of 
^O'.'l size (~70 AU) located in R5. The maser spot posi- 
tions can be fitted extremely we h by a circle o f ^^62 AU 
radius (with a precision of 0.1%. lTbrrelles et~al . 2001a), 
and their proper motions indicate uniform expansion of 
~9 km s""'^ perpendicular to the arc. Both the spatial 
distribution and the proper-motion direction of water 
masers strongly suggest that the arc is part of a spherical 
expanding structure excited by a YSO, probably loc ated 
-0'.'6 (-400 AU) south of HW2 (IGuriel et al.ll2002l l. 

W75 N is also a high mass SFR located at a distance 
of ~2 kpc (jPickel et al.l 119691 ). At scales of arcseconds 
(few thousands of AUs), there are three ultracompact 
radio continuum sources, VLA 1, VLA 2 and VLA 3 
(iHu nter et all 119941 : ITorrelles etaD 11997 ^. whic h are 
probably excited by early B -type stars (|Shcphcrd et al.l 
I2OOI I2OOI [Persi et al.ll20M ). VLA observations toward 
this region revealed two clusters of water masers, one 
associat ed with VLA 1 and the other associated with 
VLA 2 (jTorrelles et al.l 119971) . In order to study in de- 
tail not only the spatial distribution of water masers but 
also th e gas kinematics in these sources, ITorrelles et al. I 
()2003l ) carried out VLBA multi-epoch water maser obser- 
vations to measure their proper motions. They detected 
~700 maser spots in each of the three observed epochs, 
which were distributed in two groups, each one associ- 
ated with a radio continuum source, but tracing outflows 
with clearly different characteristics. In VLA 1 the maser 
spots form a linear structure of ~0'.'75 (~1500 AU) along 
the major axis of the radio continuum emission with a 
mean proper-motion value of ~19 km parallel to this 
major axis. On the other hand, in VLA 2 the maser 
spots describe a shell of ^O'.'IQ size (—320 AU) around 
the radio continuum emission. The masers move out- 
ward from the central source in multiple directions with 
a mean proper-motion value of ~28 km s~^. The de- 
gree of collimation is clearly different in the two outflows 
traced by the water masers in this SFR. 

The largest maser structures of these two SFRs are lo- 
cated in subregions R5 in Cepheus A, and VLA 1 and 
VLA 2 in W75 N. Each structure is associated with 
a YSO and could give fundamental information about 
the gas dynamics very close to the protostar. Since 
the number of maser components is very large, it is 
important to do statistical studies of the properties of 
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the spatio-kinematical distribution of these masers from 
scales smaher than 1 AU up to a few hundred of AUs. 

4. RESULTS 

4.1. Statistical properties of the spatial distribution of 
water masers 

Using VLBA data of three epochs of water maser 
observations toward the SFRs Cepheus A and W75 N 
(jTorrelles et al I l2001bL 120031 ). we compute two-point 
spatial correlation functions of maser spots in the fol- 
lowing subregions: R5 in Cepheus A, and VLA 1 and 
VLA 2 in W75 N. In each of these subregions, a few 
hundred maser spots were detected with an S/N > 10 in 
each one of the three epochs (see Table 1). The uncer- 
tainties in the values of these functions are mainly due to 
position uncertainty of maser spots, which is more signif- 
icant toward small separations (represented by the error 
bar sizes shown in Figures 1-3). We estimated the sizes of 
the error bars in the two-point spatial correlation func- 
tions in the following way: starting from the observed 
spatial distribution of maser spots, we have applied ran- 
dom displacements, following a Gaussian probability dis- 
tribution, with a dispersion equal to the positional uncer- 
tainty. We have calculated several such alternative dis- 
tributions of masers and estimated the correlation func- 
tion of each of them. In this way, we obtained a set of 
different alternative correlation functions. Then, we cal- 
culated the standard deviation of the set of values of the 
correlation function at each separation and adopted that 
deviation as the typical error of the correlation function. 

We approximate the two-point spatial correlation func- 
tions by two different power laws, one within small sepa- 
rations and the other within larger separations, because 
there is a hint of a break in the slope of the correlation 
function. An approximation with two different power 
laws has a smaller value of chi-square than a single power- 
law fit. We calculate the best value for the breakpoint 
as that for which the chi-square in the fit is minimum. 

We have adopted a minimum value of the scale range 
of 0.07 AU for Cep A R5, and of 0.2 AU for VLA 1 and 
VLA 2 in W75 N, corresponding to two times the accu- 
racy in relative positions of maser spots {2a ~ 0.1 mas) 
for each case. At separations smaller than this value, the 
correlation function is not accurate, since positional in- 
formation is dominated by random noise. The power-law 
indices and the scale ranges of the fits are listed in Table 
2. 

Two-point correlation functions for maser spots in Cep 
A R5 for the three observed epochs are shown in Figure 1. 
These functions follow a power-law dependence between 
0.07 and 0.24 AU with a weighted mean value of the index 
a = -2.79 ± 0.22 (Table 2), illustrated as a straight line 
in this figure. At a separation of 0.31±0.07 AU, there 
is a break of the power- law slope (indicated by shades). 
A different power-law dependence is found in the scale 
range of 0.5—70 AU, with a weighted mean value of the 
index, a = -1.46 ± 0.09 (Table 2). 

Figure 2 shows the two-point correlation functions for 
VLA 1 in W75 N. In VLA 1, the number of maser spots 
detected is smaller than in the other regions, introducing 
larger errors in our statistical analysis. In fact, we have 
found that the values of the correlation function at sep- 
arations smaller than 1 AU are more uncertain. Large 



error bars are due to both the relative position uncer- 
tainty of the maser spots and small number (a few tens) 
of maser spots at these separations. At separations be- 
tween 0.2 and 0.75 AU, the correlation function roughly 
follows a power with a weighted mean value of the index 
a = -3.92 ± 0.33 (Table 2). There is a hint of a break 
of the power-law slope at a separation of 0.85±0.10 AU 
(indicated by shades). Another power law is found in the 
scale range of 1.3—1500 AU, with a weighted mean value 
of the index, a = -1.47 ± 0.10. 

Figure 3 shows the two-point correlation functions for 
VLA 2 in W75 N. These functions follow a power law 
between 0.2 and 0.6 AU with a weighted mean value of 
the index a = -2.51 ± 0.17 (Table 2). There is a break 
of the power-law slope at a separation of 0.9±0.3 AU 
(indicated by shades). A different power law is found in 
the scale range of 1.5—250 AU, with a weighted mean 
value of the index, a — —1.30 ± 0.07. 

4.2. Statistical properties of the velocity distribution of 
water masers 

We have calculated the velocity correlation functions 
defined in §2.2 for water maser spots in the subregions 
R5 in Cepheus A, and VLA 1 and VLA 2 in W75 N. The 
velocity resolution was 0.21 km s^^ and 0.4 km s~^ for 
the Cepheus A and W75 N observations, respectively. 

In general, both velocity correlation functions, Vs{Ar) 
and Afs(Ar), can be approximated by a power law in 
the whole range. The power-law indices and the scale 
ranges of the fits are given in Table 3. These two veloc- 
ity correlation functions, the rms and the median of the 
difference in Doppler velocity as a function of the sep- 
aration between spots, are shown in Figure 4 for subre- 
gion Cep A R5. The obtained power-law indices for both 
functions are consistent with zero for the three observed 
epochs, with weighted mean values a = —0.03 ± 0.03 
and a = 0.00 ±0.04 for the rms and the median Doppler- 
velocity difference, respectively (Table 3). The error bars 
in the velocity correlation functions are due to both the 
uncertainty in the relative positions of the maser spots 
and the uncertainty in their velocities. These error bars 
were estimated in an analogous way to the two-point 
spatial correlation function, but including the velocity 
uncertainty, i.e., letting the velocity of the maser spots 
vary within the velocity uncertainty (the channel width 
in the observations). Then we estimate the velocity cor- 
relation function of each new distribution and compare 
the results with the original, estimating a typical error 
from the standard deviation of the velocity correlation 
function after several trials. 

Figures 5 and 6 show the velocity correlation functions 
for VLA 1 and VLA 2 in W75 N. In these subregions, the 
correlation functions show a trend of increasing Doppler- 
velocity difference with increasing separation. If these 
functions are approximated by a power law, in VLA 1 
the weighted mean value of power-law indices are a = 
0.17 ± 0.05 and a = 0.14 ± 0.05 for the rms and the 
median Doppler-velocity difference, respectively (Table 
3). In VLA 2, the weighted mean value of the power-law 
indices are a = 0.13±0.03 and a — 0.20±0.02 for the rms 
and the median Doppler-velocity difference, respectively. 

5. DISCUSSION 
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5.1. Clustering 

In order to understand the information given by the 
two-point spatial correlation functions of maser spots in 
these SFRs, we carried out a similar statistical analysis 
for some hypothetical cases with an illustrative purpose. 
For example, we calculated the two-point spatial correla- 
tion function for a cluster of ~10000 dots randomly and 
uniformly distributed in a circle with a diameter equal 
to 1.0 (in arbitrary units). The results are shown in the 
first row of Figure 7. At separations much smaller than 
1, the curve has an index a = 0.0 because the number 
of dots per unit area does not depend on dot separation. 
The curve drops rapidly at separations close to the di- 
ameter of the circle, as a result of the fact that the dots 
are confined to a finite region. 

We have also studied the case of a discrete circle with a 
diameter equal to 1.0 formed by several circular clusters, 
each one with a diameter equal to 0.15. In the case of 
three clusters, the distance between any pair of clusters 
is the same (see Figure 7, second row). The two-point 
correlation function for such a distribution shows that 
the number of dots per unit area decreases rapidly at 
separations close to the diameter of each cluster. As sep- 
aration becomes larger than 0.15, the number of dots per 
unit area increases gradually, because we start to count 
dots that belong to another cluster, having a maximum 
at a separation equals to the value of the distance be- 
tween the centers of the clusters. Increasing the number 
of clusters results in an increase of the number of max- 
ima, at separations equal to the different values of the 
distance between the centers of the clusters (see Figure 
7, third and fourth rows). 

From the hypothetical cases described above, we know 
that the two-point spatial correlation function of a ran- 
dom and uniform distribution of dots can be described by 
a power-law dependence with an index a — 0.0. There- 
fore, the power-law indices found for the maser spots in 
the SFRs Cepheus A and W75 N (see Table 2) strongly 
suggest a characteristic clustering of maser spots. This 
behavior indicates that the water maser spots are associ- 
ated in a scale-free or fractal fashion on the sky, as noted 
for other SFRs such as W49N, W3 IRS 5, and Sgr B2(M ) 
(|Gwinn Ill994at llmai et al. l[200l IStrelnitski et"all[200l . 

From Figure 7, we can estimate the characteristic size 
of the clusters of dots as the scale or separation where 
the first steep drop of the two-point correlation function 
occurs. Inspection of Figures 1, 2, and 3 reveals a break 
in the slope of the two-point correlation functions of the 
three studied regions (Cep A R5, W75 N VLA 1, and 
W75 N VLA 2) that is highlighted by shades. If we 
identify this breakpoint as the characteristic scale, we 
can estimate a typical size for clusters of maser spots 
of ~0.31±0.07 AU in Cep A R5. This value is slightly 
smaller than the value found in W3 IRS 5 ('^0.5 AU, 
llmai et al. |[200l . For VLA 1 and VLA 2 in W75 N, the 
typical size for clusters of maser spots is ~0.85±0.10 AU, 
and ~0.9±0.3 AU, respectively. These values found in 
W75 N are very close to the typical size of wat er maser 
featur es of 1 AU found in the SFR W49 N (jGwinn I 
ll994aD . We note that the typical size for clusters of wa- 
ter maser spots tends to be <1 AU in all these SFRs. 
This may be related with a scale for water maser excita- 
tion. Theoretical models of water masers based on shock 



excitation yield to a similar scale size for masing regions 
(|Elitzur et al.l[T98l Kaufman fc Neufeldi ri996'). 

Additionally, we note that in the regions studied in 
this work, the typical size of the clusters of maser spots 
is smaller in Cepheus A than in W75 N. This could re- 
flect a systematic error (for instance, in distance deter- 
mination) , or be an intrinsic property of the regions (de- 
pending on their evolutionary stage or stellar mass, for 
instance). Studies in other regions with the tools used in 
this paper would be necessary to identify whether these 
differences have a physical meaning. 

At separations larger than the characteristic scale of 
clusters, the two-point correlation functions of the mod- 
eled distribution of clusters show oscillations with local 
maxima located at separations equal to the values of the 
distance between the centers of the clusters (see Figure 
7). If we observe in detail the close-up of the two-point 
correlation functions of maser spots in the three observed 
regions (Figs. 1, 2, and 3), we then note these oscil- 
lations at separations larger than the breakpoint. The 
local maxima may correspond to the typical distances 
(-0.7 AU in Cep A R5, and ~3 AU in W75 N VLA 1 
and W75 N VLA 2) between clusters of maser spots. 

The change of the power-law index between small 
scales (a = —4.3 to —2.4) and large scales (a = —1.56 
to —1.27) mentioned in §4.1, suggests that the statistical 
properties of the clustering of maser spots are different at 
scales smaller and larger than the typical size of clusters. 

Furthermore, we have found that the value of the index 
a hardly changes from one epoch to another, at separa- 
tions larger than the characteristic scale size for clusters 
of water masers (Table 2). This behavior suggests that 
the spatial distribution of clusters of water masers does 
not change significantly on large scales during the three 
observed epochs. In fact, structures of several tens of 
AU are preserved during the time span of the observa- 
tions (two months). On the other hand, the value of the 
index a apparently changes from one epoch to another 
at separations smaller than the characteristic scale size 
for clusters (Table 2). This probably indicates a change 
in the density or spatial distribution of the maser spots 
inside each cluster. 

5.2. Velocity distribution 

IStrelnitski et al.l ()2002D investigated whether regular 
(non-turbulent) velocity fields, such as expansion and/or 
rotation can produce the power-law dependence of the 
velocity increments on spatial scales, described in the 
velocity correlation functions. Specifically, they used a 
model of 90 dots randomly and uniformly distributed in 
a thin spherical shell considering radial expansion, rota- 
tion around an axis perpendicular to the line of sight, and 
both motions. The line-of-sight velocity difference as a 
function of dot separation c an be fitted by a power law in 
a given range (see Fig. 8 of IStrelnitski et al.ll2002t) . with 
an index of about zero for the case of expansion, changing 
to approximately unity for the case of rotation. When 
both motions are present, the index has intermediate val- 
ues. Based on these results, for a spherical geometry, we 
can conclude that if organized motion (expansion and/or 
rotation) dominates over turbulence over a certain scale 
range, the value of the index is expected to be between 
0.0 and 1.0. 

In the Cep A R5 and W75 N VLA 2 regions, wa- 
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ter masers delineate remarkable arc structures, showing 
proper motions with a mean value of ~9 km s~^ and 
^28 km s~^, respectively (see §3). The spatial distri- 
bution of the masers and the direction of the proper 
motio ns are consistent with exp anding motions in a 
shell (jTorrelles et al. ll2001bL 120031 ) . On the other hand, 
in the W75 N VLA 1 region, water masers form a 
linear structure with a mean proper-motion value of 
~19 km s^^ along the direction of the elongated struc- 
ture that was interpreted as tracing a coUimated outflow 
(jjorrellcs ct al. 2003). 

In our statistical analysis of the velocity distribution of 
the water masers in Cep A R5, we find that the velocity 
correlation functions of both the rms and the median 
Doppler-velocity difference can be fitted by a power law 
with an index of about zero for the three observed epochs 
(Figure 4). This index is consistent with the case of 
pure expansion of dots randomly dis tributed in a thin 
spherical shell (jStrelnitski et al.ll2002|) . This result is in 
agreement with the VLBA proper-motion measurements 
of wa ter masers, which indic ate expansion of a spherical 
sheU (Torrelle s et"ani2001al lb[). 

On the other hand, the value of the power-law index 
found for the velocity correlation functions of maser spots 
in W75 N VLA 2 is ~0.18, suggesting another geome- 
try (different from a spherical shell) for this source. In 
fact, in W75 N VLA 2, masers delineate an incomplete 
elliptical ring on the plane of the sky; while in Cep A 
R5, masers delineate part of a circle on the plane of the 
sky (a circle in projection is produced only by a limb- 
brightened spherical shell). In addition, the masers in 
W75 N VLA 2 present a higher Doppler-velocity dis- 
persion than maser spots in Cep A R5 (25 km s~^ vs. 
2 km s-^: lTorrelles et al. I [200TH 120031) . 

For the reasons described above, we have considered a 
ring with an arbitrary orientation with respect to the ob- 
server, applying a similar statistical analysis to a model 
that simulates organized motions of masers. Our aim is 
to determine whether a particular orientation and the 
presence of expansion and/or rotation motions could ex- 
plain the trend of increasing Doppler-velocity difference 
with increasing separation of spots indicated by the value 
of the power-law index of the velocity correlation func- 
tion in W75 N VLA 2 (Fig. 6 and Table 3). 

We have used ~1000 dots randomly and uniformly dis- 
tributed in a ring with an outer radius equal to ^O'.'OB 
(~160 AU), the radius of the region traced by masers in 
W75 N VLA 2. The inner radius of the ring is chosen 
in such a way that the difference between the outer and 
inner radii is equal to ^0'.'02 (^40 AU), the spatial dis- 
persion of the masers observed within the VLA 2 shell. 
Also, the vertical thickness of the ring is assumed to have 
the same value. We have analyzed three types of orga- 
nized motions: uniform expansion, rigid-body rotation 
around an axis with an inclination angle 6 with respect 
to the line of sight, and expansion plus rotation. In our 
calculations, the magnitudes of the expansion and rota- 
tion velocities are assumed to be equal. In order to esti- 
mate the uncertainties involved in these calculations, we 
have considered that each different spatial configuration 
of the dots randomly distributed within the ring would 
result in slightly different velocity correlation functions. 
Therefore for each type of organized motion, we have an- 



alyzed different spatial configurations of dots to find the 
mean values of the Vs function and its dispersion. The 
results of the rms line-of-sight velocity difference {Vs) as 
a function of dot separation are qualitatively similar to 
the median line-of-sight velocity difference (M^); there- 
fore, we only present the Va function. 

The mean values of (14) as a function of dot separa- 
tion are shown in the top panel of Figure 8, presenting 
a comparison between the three types of motions for a 
ring observed with an inclination angle of 6* = 20° (the 
angle between the line of sight and the normal to the 
ring plane; 6 — 0° corresponds to the face-on case). The 
magnitudes of both the expansion and rotation velocities 
are equal to '^44 km s~^. These values of the inclination 
angle and velocities yield the best fit of the data. For 
the case of expansion (black dotted curve), the Vs func- 
tion is almost constant for dot separations smaller than 
^O'.'Ol (^20 AU), and increases smoothly for larger dot 
separations. This trend is not well approximated by a 
power-law dependence. For the case of rotation (black 
dashed curve), Vs increases continuously with dot sepa- 
ration. This tendency can be satisfactorily approximated 
by a power-law dependence with an index of about 1. 
For the case of expansion plus rotation (black continu- 
ous curve), Vs is similar to that of the pure expansion 
case at separations smaller than '^O'.'Ol (^20 AU), how- 
ever when including rotation, Vg increases more rapidly 
for larger dot separations. 

In the middle panel of Figure 8, we show a compari- 
son between the velocity correlation functions calculated 
for different spatial configurations of dots randomly and 
uniformly distributed in an arc, equal to one half of the 
previously studied ring with the same orientation with re- 
spect to the observer, and with equal values of expansion 
and rotation velocities. The selected section of the ring is 
predominantly redshifted. The values of Vs as a function 
of dot separation are similar to those for a complete ring 
for separations <0'.'05 (<100 AU); at larger separations, 
Vs decreases slightly with respect to the complete ring 
case, because in the case of a half ring the maximum ve- 
locity difference is smaller. Since maser spots in W75 N 
VLA 2 do not trace a continuous or complete structure 
on the plane of the sky, we consider that this scenario 
might help us to understand the decrease of Vs observed 
in the velocity correlation functions on a scale larger than 
-100 AU (see Figure 6). 

In fact, the bottom panel of Figure 8 shows the rms val- 
ues of Doppler-velocity difference {Vs) calculated for the 
masers in W75 N VLA 2 observed in the second epoch 
(when more maser spots were detected), together with 
the modeled functions that best fit these data (shown 
in the middle panel). Between 1 and 100 AU (verti- 
cal dotted lines), there is a general agreement between 
the functions calculated for the modeled spots and the 
observed masers. For the modeled functions, consider- 
ing pure expansion or expansion plus rotation produce 
similar results in this scale range. However, at separa- 
tions larger than 100 AU, the Vs function of the modeled 
spots decreases more rapidly when rotation motions are 
included, which is qualitatively similar to the Vs of the 
observed masers at those separations. At separations 
<1 AU, a larger dispersion is present in the velocity cor- 
relation functions for the modeled spots (shown by the 
grey curves), due to a lack of clustering, and thus poor 
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statistics. This occurs because even though the total 
number of modeled spots is large, the number of pairs 
with small separation is small because modeled spots do 
not form clusters (they are uniformly and randomly dis- 
tributed). 

From our results, we thus conclude that the trend de- 
scribed by the rms Doppler-velocity difference of maser 
spot pairs in W75 N VLA 2 as a function of separation 
can be explained by the presence of organized motions, 
expansion, and rotation (both velocities with the same 
magnitude '-~^44 km s~^) within an arc of a circle ob- 
served with small inclination angle, 9 ~20°. Higher val- 
ues of the angle 9 will produce a steeper gradient of the 
velocity difference toward larger separations. Moreover, 
the corresponding projected ellipse on the plane of the 
sky would have its minor axis much smaller than its ma- 
jor axis, which is not observed. The expansion velocity 
required to fit the data is about the sa me order of mag- 
nitud e as the measured proper motions (jTorrelles et al. 1 
[200l . 

The source W75 N VLA 1 appears to be different 
from VLA 2 since proper-motion measurements of wa- 
ter masers indicate a collimated outflow in this region 
(jTorrelles et al. II2003D . In fact, the velocity correlation 
functions of maser spots present a larger dispersion from 
the fit. Despite the dispersion, the values estimated for 
the indices (similar to those of VLA 2) suggest the pres- 
ence of an organized component of motion with a pref- 
erential direction. 

6. CONCLUSIONS 

In this paper, we present a statistical analysis of the 
spatial and velocity distribution of water masers in the 
SFRs Cepheus A and W75 N, using the data of the 
VLB A maser observations carried out bv lTorrelles et al71 
(j2001bl 12001 . Our conclusions are as follows: 

1. We have found a characteristic scale size for clus- 
ters of water maser spots <1 AU, indicated by the 
break in the slope of the two-point spatial corre- 
lation function. Specifically, ^0.31±0.07 AU in 
Cep A R5, -0.85±0.10 in'W75 N VLA 1, and 
~0.9±0.3 AU in W75 N VLA2. These values are 
close to the typical sizes of the water maser features 
found in other SF Rs, e.g., -0.5 AU in W 3 IRS 5 
(llmai et al. | [200^ and -1 AU in W49 N (jGwinn I 
ll994aD . Probablv. these results may indicate that 
the scale for water maser excitation tends to be 
<1 AU, as it was pointed ou t by models of wa- 
ter masers excited by sho cks (jElitzur et al.l 119891 : 
Kaufman fc Neufeldlll996l) . 

2. Two-point spatial correlation functions follow 
power-law dependences, indicating self-similar spa- 
tial distributions of water masers. The power-law 
indices found for separations smaller than the char- 
acteristic scale size of clusters (a = —4.3 to —2.4) 
are steeper in Cepheus A and W75 N than the 
value found in W3 IRS 5 (a = -2.09. llmai et al. I 
I2002D . This could be due to a different spatial dis- 
tribution of spots inside the clusters of each region. 
The power-law indices found for separations larger 
than the characteristic scale size (a = —1.56 to 
— 1.27) are similar to the value found in W49 N 



(a = -1.33. [Gwhm1ll994aj) . This suggests simi- 
lar statistical properties for the spatial distribution 
of masers at scales larger than 1 AU in the SFRs 
Cepheus A, W75 N, and W49 N. 

3. Velocity correlation functions follow power-law de- 
pendences. In Cep A R5, the power-law index is 
about zero, consistent with pure expansion of dots 
randomly distributed in a thin spherical shell. In 
W75 N VLA 2, the value for the power-law in- 
dex can be explained by the presence of expan- 
sion and rotation within an arc of a circle observed 
with a small inclination angle. These results are 
in agreement with proper-motion ob ser vatio ns of 
water masers (jTorrelles et al. Il2001al lbl. I2003D . In 
W75 N VLA 1, the value for the power-law index 
seems to be consistent with organized motions but 
further studies are needed. In the SFRs Cepheus A 
and W75 N, the values estimated for the power-law 
indices suggest that water masers are tracing an or- 
ganized component of the velocity field instead of 
a turbulent component, probably due to a preva- 
lence of organized motions over turbulent motions 
at spatial scales from a few hundred of AUs down 
to less than 1 AU. 

4. Based on the value of the power-law index followed 
by the velocity differences of maser spot pairs as 
a function of their separation, it may be possible 
to find evidence for expansion and / or rotation mo- 
tions by analyzing water maser data from a single 
observed epoch. 
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Fig. 1. — Left panels: Two-point correlation functions for water mascr spots in subrcgion Ccp A R5 for cacii of tiic three observed epochs. 
The straight lines show power-law fits in two different scale ranges, 0.07—0.24 AU and 0.5—70 AU (see Table 2). The shaded region shows 
the scale range within which the power-law slope changes. The positional accuracy (2(t) of the spots is 0.07 AU. The error bars arc mainly 
due to position uncertainty of maser spots. Right panels: Close-up of the correlation function for each of the three observed epochs. We 
recomputed this function using a smaller value of dr to see in detail the break of the power-law slope. The lines show the same power-law 
fits as in the left panels. 
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Fig. 2. — Left panels: Two-point correlation functions for water maser spots in subrcgion W75 N VLA 1 for each of the three observed 
epochs. The straight lines show power-law fits in two different scale ranges, 0.2—0.75 AU and 1.3—1500 AU (see Table 2). The shaded 
region shows tlic scale range within which the power-law slope changes. The positional accuracy (2<t) of the spots is 0.2 AU. Right panels: 
Close-up of the correlation function for each of the three observed epochs. We recomputed this function using a smaller value of dr to see 
in detail the break of the power-law slope. The lines show the same power-law fits as in the left panels. 
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Fig. 3. — Left panels: Two-point correlation functions for water maser spots in subregion W75 N VLA 2 for each of the three observed 
epochs. The straight lines show power-law fits in two different scale ranges, 0.2—0.6 AU and 1.5—250 AU (see Table 2). The shaded 
region shows tlic scale range within which the power-law slope changes. The positional accuracy (2<t) of the spots is 0.2 AU. Right panels: 
Close-up of the correlation function for each of the three observed epochs. We recomputed this function using a smaller value of dr to see 
in detail the break of the power-law slope. The lines show the same power-law fits as in the left panels. 
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Fig. 4. — Velocity correlation functions for water maser spots in Cep A R5 for each of the three observed epochs. Left panels show the 
rms values of Doppler- velocity differences (Vs) between maser spot pairs as a function of their separation on the plane of the sky. Right 
panels show the medians of the Doppler- velocity differences (Ms) between maser spot pairs as a function of their separation on the plane 
of the sky. The straight lines show power-law fits (Table 3). The positional accuracy (2<t) of the spots is 0.07 AU. Error bars are due to 
position and velocity uncertainties of maser spots. 
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Fig. 5. Velocity correlation functions for water maser spots in W75 N VLA 1 for each of the three observed epochs. Left panels show 
the rms values of Dopplcr- velocity differences (Vs) between maser spot pairs as a function of their separation on the plane of the sky. Right 
panels show the incdiaiis of tiic Dopplcr- velocity differences (Ms) between maser spot pairs as a function of their separation on the plane 
of the sky. The straight lines show power-law fits (Table 3). The positional accuracy {2a) of the spots is 0.2 AU. 
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Fig. 6. Velocity correlation functions for water maser spots in W75 N VLA 2 for each of the three observed epochs. Left panels show 
the rms values of Dopplcr- velocity differences (Vs) between maser spot pairs as a function of their separation on the plane of the sky. Right 
panels show the incdiaiis of the Dopplcr- velocity differences (Ms) between maser spot pairs as a function of their separation on the plane 
of the sky. The straight lines show power-law fits (Table 3). The positional accuracy {2a) of the spots is 0.2 AU. 
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Fig. 7. — Left panels show different distributions of dots and right panels show their corresponding two-point correlation functions. From 
top to bottom, first row represents a cluster of ~10000 dots randomly and uniformly distributed in a circle with a diameter equal to 1.0. 
Second, third, and fourth rows represent groups of S, 5, and 11 circular clusters of dots with a diameter of 0.15, respectively. Dots are 
randomly distributed within each cluster. Clusters are equally spaced along a circle with diameter equal to 1.0. In right panels, the first 
dotted line indicates a value of the separation equal to the diameter of a cluster (0.15), and the following dotted lines indicate the distances 
between the centers of each pair of clusters. 
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Fig. 8. — Top: Velocity correlation functions for ~1000 dots randomly and uniformly distributed in a ring with an inner radius of 120 AU 
and an outer radius of 160 AU, observed with an inclination angle 8 = 20°. The black dotted, dashed, and continuous curves show the 
mean values of Vs as a function of separation for models of organized velocity fields: expansion, rotation, and expansion plus rotation. 
The grey curves show the 1<t error associated witii these simulations. The magnitudes of both the expansion and rotation velocities are 
equal to 44 km s^^. The correlation functions are not calculated at separations smaller than 1 AU, because of the increased dispersion (sec 
details in §5.2). Middle: Same as top, but only for one half of the ring using the same kinematical parameters and orientation. Bottom: 
Comparison between the rms values of Doppler- velocity differences (Vs) for maser spots in W75 N VLA 2 detected in the second epoch 
(filled circles) and the velocity correlation functions for pure expansion and expansion plus rotation, that correspond to half a ring, shown 
in the middle panel. 
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TABLE 1 

Water maser spots detected in the subregions Cep A R5, W75 N VLA 1, and W75 N VLA 2. 





Cep A R5 




W75 N VLA 1 


W75 N VLA 2 


Epoch 




^spot* 


Epoch 




Epocli "spot" 


1996 Feb 11 




283 


1999 Apr 2 


126 


1999 Apr 2 521 


1996 Mar 10 




170 


1999 May 7 


174 


1999 May 7 791 


1996 Apr 13 




234 


1999 Jun 4 


134 


1999 Jun 4 613 



* Number of maser spots detected with an S/N > 10. 



TABLE 2 

Best-fitting power laws to two-point correlation functions. 





Cep A R5 




W75 N VLA 1 


W75 N VLA 2 




index a " 


index a " 


Index a " 


Index a " 


Index a " 


Index a 




small scales 


large scales 


small scales 


large scales 


small scales 


large scales 


Epoch 


(0.07-0.24 AU) 


(0.5-70 AU) 


(0.2-0.75 AU) 


(1.3-1500 AU) 


(0.2-0.6 AU) 


(1.5-250 AU) 


1st 


-2.8±0.6 


-1.48±0.15 


-2.9±0.9 


-1.47±0.16 


-2.6±0.4 


-1.27±0.13 


2nd 


-2.4±0.4 


-1.56±0.25 


-4.3±0.4 


-1.49±0.18 


-2.44±0.20 


-1.37±0.12 


3rd 


-3.0±0.3 


-1.41±0.14 


-3.2±0.8 


-1.45±0.16 


-2.8±0.5 


-1.28±0.10 



Note. — Uncertainties are 2(t. 

Two-point correlation functions for maser spots can be fitted by: ns(Ar) = noAr" 



TABLE 3 

Best-fitting power laws to velocity correlation functions. 



Epoch 


Cep A R5 ^ 




W75 N VLA 1 ^ 


W75 N VLA 2 " 


Index a 


Index a ° 


Index a 


Index a ° 


Index a 


Index a ° 


1st 


-0.04±0.06 


-0.01it0.06 


0.13±0.07 


0.15±0.08 


0.14±0.06 


0.19±0.05 


2nd 


-0.02±0.05 


O.OlitO.06 


0.14±0.09 


0.14±0.09 


0.12±0.04 


0.16±0.04 


3rd 


-0.03±0.06 


0.02±0.08 


0.31±0.11 


0.14±0.08 


0.15±0.05 


0.23±0.04 



Note. — Uncertainties are 2a. 
a Scale range: 0.07-70 AU. 
^ Scale range: 0.2-1500 AU. 

Scale range: 0.2-320 AU. 

Index of the power-law fit to the rms velocity difference between maser spot pairs, Vs(Ar) oc Ar". 
" Index of the power-law fit to the median velocity difference between maser spot pairs, Ms(Ar) oc Ar". 



